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— ‘Doppler  observations  of  Navy  Navigation  Satellites  have  yielded  precisions 
in  the  determination  of  rates  of  change  of  station  position  of  2 cm/yr,  based 
on  data  sampled  over  a 14-yr  period,  and  3 cm/yr  for  nearly  continuous  data  over 
a 6-yr  period.  The  accuracy  of  the  determination  is  better  than  10  cm/yr  in 
horizontal  position  and  15  cm/yr  in  height  for  either  data  span.  A number  of 
changes  in  data  acquisition  and  data  reduction  are  proposed  which  will  improve 
the  accuracy  of  the  results  and  lead  to  the  determination  of  plate  tectonic 
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Doppler  observations  of  Navy  Navigation  Satellites  have  been  used 
to  determine  the  positions  of  navigation  beacons,  interdatum  relations 
and  geocentric  positions  of  isolated  sites,  and  major  datums  since  1964. 
The  calculations  have  also  provided  the  motion  of  the  earth's  spin  axis 
with  respect  to  the  crust  since  1969.  Periodically,  coordinates  of  the 
worldwide  base  stations  used  to  compute  the  satellite  orbits  are  deter- 
mined from  those  orbits  (a  circular  kind  of  computation)  to  determine 
the  stability  of  the  coordinate  system  and  identify  inconsistencies  in 
the  computational  procedure.  This  report  examines  the  results  of  these 
computations  to  evaluate  the  possibility  of  using  the  data  to  determine 
plate  tectonic  motion. 

The  precisions  achieved  to  date  in  the  calculation  of  positions  of 
Doppler  observing  stations  are  sufficient  for  determining  plate  tectonic 
motion,  but  inaccuracies  in  the  collection  and  computational  procedures 
render  the  results  unreliable.  It  is  suggested  that  continued  data  col- 
lection under  controlled  conditions  and  recomputation  of  station  posi- 
tions with  data  collected  to  date  could  provide  useful  results  on  plate 
tectonic  motion. 

Satellite  ephemerides  used  in  this  study  and  observational  data 
have  been  provided  by  Caroline  Leroy  of  the  Defense  Mapping  Agency 
Topographic  Center  since  1973.  Computation  of  station  positions  was 
made  by  Teck  Judd,  while  station  coordinate  changes  were  computed  by 
Natalie  Krohn,  both  of  the  Naval  Surface  Weapons  Center. 
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INTRODUCTION 


The  Navy  Navigation  Satellite  Program  was  established  in  the  late 
1950s  to  provide  the  real-time  position  of  ships  to  100-m  accuracy.  As 
this  goal  was  being  approached,  it  became  evident  that  post  processing 
of  data  could  be  used  for  geodetic  applications.  A project  to  develop 
this  capability  was  established  in  1962  with  a goal  of  10-m  accuracy. 
Positions  of  navigation  beacons,  isolated  tracking  stations,  datum 
origins,  and  interdatum  ties  were  found  to  50-m  accuracy  in  1963,  im- 
proving to  about  1-m  accuracy  in  1972. 1 The  improvement  in  accuracy  was 
achieved  primarily  by  refinements  in  the  determination  of  the  effects 
of  the  earth's  gravity  field  on  satellite  orbits,  except  for  the  change 
in  1969  to  the  real-time  determination  of  the  motion  of  the  earth's  spin 
axis  with  respect  to  the  crust.2  Since  that  time,  changes  in  the  com- 
putational model  have  primarily  centered  about  the  development  of  con- 
sistent station  coordinates  for  the  observing  stations. 


In  1968,  a review  of  results  to  that  time  indicated  that  20  yr 
of  intermittent  or  continous  determinations  of  station  positions  would 
yield  continental  drift  to  an  accuracy  of  15  cm/yr  or  5 cm/yr,  respec- 
tively. ^ Since  only  preprints  of  the  test  of  that  report  were  distrib- 
uted, it  is  reproduced  as  Appendix  A for  historical  perspective  and  to 
document  the  procedure  for  establishing  the  effective  position  of  the 
antenna  array.  A subsequent  review  showed  that  projection  to  be  pessi- 
mistic due  to  refinements  incorporated  in  the  calculations  in  the  interim 
period.4  A standard  error  of  10  cm/yr  was  found  for  a period  of  2.5  yr, 
and  an  accuracy  of  2 cm/yr  was  predicted  for  continuous  observations 
over  a period  "longer  than  5 years"  (7  yr  for  the  data  given).  Contin- 
uous computations  of  positions  have  now  been  completed  at  nine  locations. 
The  purpose  of  this  report  is  to  reexamine  the  potential  of  Doppler 
observations  of  Navy  Navigation  Satellites  to  determine  plate  tectonic 
motion  on  the  basis  of  these  computations  and  others  that  were  made  to 
monitor  the  stability  of  the  coordinate  system  and  detect  computa- 
tional inconsistencies. 


OBSERVING  NETWORK 

Results  of  computing  the  positions  of  22  of  the  sites  shown  in 
Figure  1 (Appendix  B)  are  discussed  in  this  report.  The  break  in  oper- 
ation of  the  stations  shown  in  the  figure  for  the  1970-1971  time  period 
was  caused  by  the  change  in  the  equipment  to  count  Doppler  in  contig- 
uous intervals  of  about  20  sec  in  length  rather  than  a 1-sec  interval 
starting  each  4 sec.  This  change  required  a change  in  data  representa- 
tion from  frequency  to  range  difference.  (The  break  in  operations  was 
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shorter  than  that  shown  in  the  sketch. ) The  22  sites  are  all  those 
shown  operating  in  the  1972.0-1978.0  time  period,  except  Grasse  and  Wake 
operated  for  only  a brief  period.  The  figure  shows  that  antenna  changes 
were  made  at  six  of  these  sites  during  that  time  period  and  that  station 
operations  were  not  conducted  during  a substantial  protion  of  the  period 
for  six  additional  sites.  In  addition,  antennas  are  switched  intermit- 
tently at  three  of  the  sites  (310,  320,  and  330).  Therefore,  a contin- 
uous record  of  data  is  available  from  only  eight  sites  for  which  no 
record  of  an  antenna  change  has  been  recorded  during  the  6-yr  period. 
However,  survey  records  between  old  and  new  antennas  at  stations  116, 
310,  320,  and  330  appear  to  be  quite  reliable,  and  the  difference  be- 
tween the  station  13  and  27  antennas  was  based  on  5 wk  of  simultaneous 
satellite  observations  at  the  two  sites.  Therefore,  it  is  believed  that 
the  continuous  record  of  observations  may  be  reliable  for  at  least  14 
sites  during  this  6-yr  period. 


COMPUTATIONAL  PROCEDURES 

Satellite  orbit  computations  were  based  upon  a least-squares  solu- 
tion for  six  constants  of  integration  of  the  orbit  of  one  of  the  Navy 
Navigation  Satellites,  an  atmospheric  drag  scaling  factor,  and  the 
components  of  pole  position  using  2 days  of  observations  from  the  ob- 
serving network.  (In  a few  of  the  computations  at  the  beginning  of  the 
1964-1978  time  period,  longer  time  spans  of  data  were  used,  in  which 
case  additional  drag  parameters  were  included.)  The  parameter  set  also 
included  a refraction  scaling  parameter  and  oscillator  frequency  offset 
for  each  pass  and  the  coordinates  of  new  observing  stations.  The  gravity 
field  was  defined  by  the  NWL-9B  geopotential  coefficients  until  2 Jan- 
uary 1973,  at  which  time  the  NWL-10E  coefficients  were  used.  The  coor- 
dinates of  base  stations  were  nominally  defined  as  the  NWL-9D  set  until 
October  1977  when  the  NWL-9Z2  set  was  introduced.  However,  each  of 
these  coordinate  sets  was  extended  to  include  values  for  new  stations, 
as  they  became  available,  and  occasionally  to  correct  previous  additions. 
Computations  for  observations  made  before  1972  were  generally  done  in 
the  1970-1972  time  period  in  connection  with  the  revision  of  the  Depart- 
ment of  Defense  World  Geodetic  System  1972.  The  remaining  computations 
produced  the  precise  ephemerides  used  by  the  Defense  Mapping  Agency  and 
other  agencies  for  percise  positioning.  The  same  observations  used  to 
determine  the  satellite  ephermeris  were  then  used  to  determine  the 
position  of  the  ground  station  with  respect  to  the  ephemeris.  Five  days 
of  observations  are  used  to  determine  the  position  of  the  station.  Also 
included  in  the  solution  was  a scale  factor  for  tropospheric  refraction 
and  the  oscillator  frequency  offset  for  observations  made  after  October 
1975  and  for  results  given  for  1972.  In  the  computations  made  prior  to 
October  1975,  the  tropospheric  refraction  was  held  fixed.  The  1973 
results  in  this  report  were  computed  after  the  computer  program  was 
changed  to  include  the  refraction  parameter. 
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In  1977,  errors  were  found  in  the  calculation  of  the  effect  of  the 
travel  time  of  the  signal  from  the  satellite  to  the  receivers.  The 
error  did  not  exist  in  the  programs  used  to  compute  the  satellite  ephem- 
eris  but  did  exist  in  the  two  programs  used  at  the  Naval  Surface  Weapons 
Center  (NSWC)  to  compute  ground  station  position  using  an  earth-fixed 
ephemeris.  To  compensate  for  the  error,  the  following  corrections  were 
added  to  the  computed  longitudes  of  stations : 

Doppler  stations:  +'.'097 
Tranet  stations*:  -'.'009 
Geoceivers**:  +'.'030 

The  computer  program  used  to  process  geoceiver  data  was  corrected  11 
May  1977  at  NSWC.  The  Geodetic  Survey  Squadron  of  the  Defense  Mapping 
Agency  Topographic  Center,  which  uses  a program  based  on  the  same  com- 
puter code,  also  tested  the  effect  of  the  program  change  and  found  the 
revised  longitudes  to  be  0. 79+. 06  m larger  for  station  latitudes  rang- 
ing from  6°S  to  65°N,  which  is  a slightly  smaller  correction.5 

A linear  least-squares  fit  was  made  to  the  5-day  solutions  for  each 
component  of  position  of  each  station  for  the  period  1972.0-1978.0  and 
for  other  cases  of  interest.  The  5-day  solutions  were  weighted  in 
accordance  with  the  number  of  passes  in  the  solution.  Solutions  which 
were  based  on  fewer  than  four  passes  were  excluded  from  the  least-squares 
fit.  Three  iterations  of  a filter  were  made  to  reject  solutions  that 
departed  from  the  linear  fit  by  more  than  2.5  times  the  standard  devia- 
tion of  the  residuals  of  fit. 


RESULTS 

Table  1 (Appendix  C)  shows  the  changes  in  coordinates  of  12  of  the 
14  stations  for  which  continuous  data  is  available  over  almost  all  of 
the  6-yr  period  1972.0-1978.0.  Results  for  Maine  and  Minnesota  are  ex- 
cluded from  the  table.  Large  rates  obtained  for  one  or  more  components 
of  position  of  these  two  stations  may  be  due  to  unusually  large  errors 
in  the  nominal  coordinates  of  the  station  used  in  the  solution  or  to 
improper  handling  of  results  from  different  antennas  used  at  the  stations. 
The  standard  errors  of  the  rates  are  close  to  those  predicted  from 
earlier  tests  with  shorter  data  spans,  particularly  since  only  part  of 


* Doppler  tracking  system  (TRANET)  is  a network  of  receiving  stations 
and  other  facilities  used  to  acquire  and  process  Doppler  data  from 
near-earth  satellites. 

**  Unless  processed  in  the  computer  program,  used  to  process  TRANET  data. 
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the  data  in  1972  was  processed  for  the  stations  and  no  data  for  1973 
was  processed  for  three  of  the  stations  listed  in  Table  1 and  for  the 
two  stations  excluded  from  the  table.  The  rates  of  change  of  latitude 
are  reasonable  in  most  cases.  Even  the  rate  of  change  of  Japan  and  the 
Phillipines  with  respect  to  Australia  is  in  the  expected  direction.  The 
rate  of  change  of  the  North  American  station  is  consistent  in  direction 
and  magnitude  with  the  drift  of  the  mean  pole  observed  classically. 
However,  the  longitude  and  height  rates  are  unreasonable  for  a large 
percentage  of  the  stations,  as  they  were  for  shorter  data  spans. 

Since  the  anomalous  rates  in  longitude  and  height  are  believed  to 
be  due  primarily  to  the  use  of  a fixed  tropospheric  model  in  the  calcu- 
lation of  refraction  effects  for  data  spans  prior  to  1972  and  in  the 
1973-1975  time  period,  the  rates  were  recomputed  using  positions  computed 
for  1972,  1976,  and  1977  data.  Results  of  these  calculations  are  given 
in  the  third  column  of  rates  in  Table  2.  Results  are  also  given  for  all 
stations  based  on  the  data  spans  1973.0-1975.5  in  the  first  column  as 
reported  by  Anderle1*,  based  on  1972-1978  data  span  in  the  second  column 
(from  which  Table  1 was  abstracted) , and  based  on  the  1964-1978  data 
span  in  the  fourth  column  of  rates.  The  increase  in  data  span  from  2.5 
to  6 yr  (columns  1 and  2)  generally  resulted  in  reduced  rates  in  latitude 
and  height  rates,  although  there  are  exceptional  cases  such  as  the 
latitude  rates  for  the  Alaskan  station  and  the  height  rates  for  the 
European  and  Maryland  stations.  The  longitude  rates  increased  about  as 
often  as  they  decreased  with  the  increase  in  span.  The  omission  of  the 
1973-1975  data  (coluirns  2 vs  3)  did  not,  in  general,  significantly  alter 
the  magnitudes  of  the  rates  compared  with  the  standard  errors.  The 
increase  in  data  span  from  6 yr  to  include  the  scattered  data  back  to 
1964  (colums  2 and  4)  generally  resulted  in  a reduction  of  uh’  rates. 

Only  a few  longitude  rates  for  the  1964-1978  span  exceed  10  cm/yr,  and 
this  is  considerably  better  than  was  expected  in  1968,  as  cited  in  the 
introduction.  However,  the  use  of  the  long  data  span  introduces  many 
additional  changes  in  antenna  locations,  for  which  documentation  becomes 
poorer  for  earlier  time  periods  when  accuracy  requirements  were  less 
stringent  and  which  affect  the  result  for  individual  stations.  But  more 
importantly,  the  variations  in  length  of  occupation  of  individual  stations 
make  comparison  of  rates  of  little  value,  since  either  geophysical  effects 
on  station  locations  or  artificial  effects  of  anomalies  in  computational 
procedures  will  yield  different  relative  rates  for  stations  that  observe 
in  different  time  periods.  Table  3 is  extracted  from  Table  2 to  summarize 
the  results  for  those  stations  that  were  occupied  a substantial  portion 
of  the  14-yr  period  from  1964-1978.  But  the  table  shows  wide  variation 
in  the  number  of  solutions  either  available  or  processed  for  these 
stations.  The  Maryland  station  was  closed  before  the  end  of  the  period, 
and  the  Alaska  and  Greenland  antennas  were  moved  to  unsurveyed  locations 
before  the  end  of  the  period,  while  1973  data  was  not  processed  for  four 
other  stations.  There  is  uncertainty  about  changes  in  antenna  locations 
at  several  of  the  sites. 
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A puzzling  result  is  the  negative  bias  in  the  longitude  rates  for 
the  1972-1978  data  spans  (Tables  1 or  2a) . If  the  computations  were 
done  properly,  such  a bias  could  not  occur  because  the  computed  satel- 
lite orbit  is  fitted  to  the  observations  with  the  station  longitudes 
held  fixed,  and  then  the  station  longitudes  are  computed  from  that  orbit. 
The  only  explanation  for  the  result  that  comes  to  mind  is  that  the  error 
in  the  positioning  programs  is  gradually  changing  the  coordinate  system 
as  new  sites  are  added  to  the  network  and  as  positions  are  computed  for 
new  antenna  locations  for  old  sites.  But  despite  this  and  other  problems 
cited  above,  the  computed  positions  of  stations  based  on  5 days  of 
observation  of  one  satellite  are  still  obtained  with  a consistency  of 
103  cm,  76  cm,  and  116  cm  in  longitude,  latitude,  and  height,  respec- 
tively, as  shown  in  Table  4 for  the  6-yr  period  1972-1978.  The  consis- 
tency over  the  14-yr  period  1964-1978  is  similar.  Furthermore,  the 
precision  and  accuracy  being  obtained  for  geophysical  applications  are 
considerably  better  than  forecast  previously.  Since  it  is  unlikely 
that  measured  rates  of  change  of  station  coordinates  would  be  too  small 
due  to  observational  or  computational  error,  the  actual  rates  provide 
a worst-case  estimate  of  the  accuracy.  Table  1 shows  that  6 yr  of  nearly 
continuous  data  has  yielded  precisions  of  3.5,  3.0,  and  3.6  cm/yr  in  rate 
of  change  of  station  longitude,  latitude,  and  height,  respectively,  and 
accuracies  of  better  than  12.5,  8.4,  and  18.8  cm/yr,  respectively. 

Table  3 shows  that  14  yr  of  intermittent  data  yields  precisions  of  1.9, 
1.7,  and  1.6  cm/yr  in  rate  of  change  of  station  longitude,  latitude, 
and  height,  respectively,  and  accuracies  of  better  than  11.0,  8.1,  and 
15.1  cm/yr,  respectively.  The  precisions  are  optimistic,  particularly 
for  very  short  data  spans,  because  they  neglect  nonlinear  effects  or 
errors  that  are  correlated  with  time.  Correlation  of  the  station 
coordinate  changes  may  be  seen  by  inspecting  the  individual  solutions 
that  are  given  in  Figure  2 for  each  site.  On  the  other  hand,  both 
precision  and  accuracy,  particularly  in  height  and  longitude,  can  be 
improved  by  improved  processing  techniques  and  by  other  methods  described 
below. 


USEFULNESS  OF  THE  DOPPLER  TECHNIQUE 

With  the  potential  of  the  Very  Long  Baseline  Interferometer  (VLBI) , 
Lunar  laser  and  Lageos  laser  techniques  to  determine  plate  tectonic  motion 
to  high  accuracy  after  a brief  period  at  a site,  it  is  reasonable  to 
question  whether  any  special  effort  should  be  pursued  to  accomplish  the 
same  task  using  Doppler  data  when  5 yr  or  more  of  nearly  continuous 
Doppler  data  would  be  required.  A number  of  factors  favor  the  duplication 
of  effort.  Of  course,  the  usual  desire  exists  for  independent  verifi- 
cation of  a difficult  measurement.  Second,  the  relatively  low  cost  of 
the  equipment  and  relatively  lower  skill  level  required  in  operation 
could  encourage  more  international  participation,  which  would  make 
results  available  from  regions  which  might  otherwise  be  inaccesible  or 
expensive  to  occupy.  Third,  weather  in  many  regions  is  not  suitable  for 
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laser  operations.  Finally,  Doppler  data,  little  of  which  has  been  pro- 
cessed for  this  purpose,  already  exists  for  sites  on  many  of  the  plates. 
Figure  3 shows  the  sites  discussed  in  this  report,  the  rates  and  stan- 
dard errors  of  rates  of  station  motion  calculated  from  the  1972-1978 
span  of  data,  * id  the  zones  of  upwelling  and  downwelling  sketched  from 
various  reports.  In  addition  to  these  sites,  geoceiver  data  exists  in 
the  Defense  Mapping  Agency  library,  and  Doppler  data  has  been  taken  by 
agencies  in  many  other  countires  at  a variety  of  locations.  Therefore, 
results  could  be  obtained  earlier  for  numerous  geographic  locations 
using  Doppler  data  than  using  the  more  precise  data. 


POTENTIAL  IMPROVEMENTS  IN  THE  DOPPLER  RESULT 

Many  steps  can,  and  some  will,  be  taken  to  improve  the  precision 
and  accuracy  of  determining  plate  tectonic  motion  from  Doppler 
observations  of  Navy  navigation  satellites.  These  will  be  discussed 
in  the  order  listed  in  Table  5 from  the  point  of  view  of  their  signifi- 
cance with  respect  to  improving  current  results,  the  cost  of  implementing 
and  operating  under  the  revision,  whether  the  change  will  be  made  in  the 
normal  course  of  events  or  need  an  administrative  decision  to  bear  the 
additional  cost,  and  whether  technical  studies  are  required  to  assess 
impact  before  the  decision  is  made.  Costs  of  changes  in  the  orbit  or 
positioning  techniques  are  estimated  as  low,  moderate,  or  high  compared 
with  the  normal  cost  of  the  respective  computations,  while  the  other 
costs  are  rated  relative  to  total  computational  costs. 

Most  of  the  base  stations  are  equipped  with  rubidium  oscillators, 
but  most  of  the  geoceivers  and  other  equipment  are  equipped  with  crystal 
oscillators.  Figure  4 shows  the  oscillator  stability  as  a function  of 
time  span  of  measurement  for  geoceiver  oscillator  specifications,  JMR* 
specifications,  measured  values  on  a cesium  oscillator  at  NSWC,  and  a 
rubidium  oscillator  at  the  Jet  Propulsion  Laboratory.  (The  stability  of 
the  geoceiver  and  JMR  oscillators  was  extrapolated  from  the  point  values 
shown  parallel  to  Allen  variance  curves  for  other  oscillators. ) Reading 
the  slanted  grid  lines  labeled  at  the  top  of  the  chart  for  averaging  times 
of  5 to  10  min,  it  is  evident  that  errors  of  20  to  40  cm  could  build  up 
in  range  difference  measurements  made  with  crystal  oscillators.  But 
this  is  the  critical  time  interval  for  observations  of  satellites  at 
altitudes  of  about  1000  km  because  the  integrated  deviation  of  the  actual 
frequency  over  this  time  period  from  the  frequency  solution  for  the  20- 
min  pass  will  directly  alias  into  the  solution  for  the  coordinates  of  the 
station  with  an  amplification  factor  of  greater  than  5.  (A  random  error 
of  observation  over  30  min  of  10  cm  corresponds  to  about  a 50-cm  standard 
error  in  station  position.  It  is  believed  that  a systematic  error  may  be 
more  serious.)  The  oscillators  are  probably  better  than  the  specifica- 
tions given,  but  the  quartz  oscillators  appear  to  be  a significant  source 
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6 


r 


of  error,  while  the  rubidium  oscillator  would  produce  errors  smaller  than 
those  due  to  jitter  in  the  phase  measurement  (11  cm  in  the  geoceiver 
specifications) . While  the  error  due  to  the  quartz  oscillator  is  still 
not  large  compared  with  the  orbit  uncertainty  under  good  geometry  con- 
ditions, a large  percentage  of  passes  occur  at  the  lower  elevation  angles 
where  the  instrument  error  dominates  the  orbit  error.  Improvement  in 
receiver  accuracy  would  mean  fewer  satellite  passes  would  have  to  be  ob- 
served to  average  out  instrument  error.  A study  is  being  conducted  to 
quantify  this  conclusion. 

Availability  of  data  from  additional  stations  by  other  agencies 
and  countries  would  be  a considerable  aid  in  isolating  errors  and 
densifying  results.  But  the  most  important  factor  is  maintaining  a 
fixed  base  network  without  changes  in  antenna  locations.  Since  there 
is  no  current  program  for  Doppler  determination  of  plate  tectonic 
motion,  station  deployment  is  determined  solely  on  the  basis  of  the 
requirements  of  the  mission  being  conducted.  Therefore,  any  program 
established  would  have  to  bear  the  cost  of  continued  operation  of  base 
stations  not  required  for  other  purposes.  A large  number  of  base 
stations  are  desirable  to  average  out  the  effects  of  motion  of  individual 
plates  on  the  computed  orbits. 

The  principal  error  in  the  satellite  orbit  is  dre  to  uncertainties 
in  the  gravity  field.  Refinements  in  the  field  are  being  made  with  the 
use  of  GEOS- 3 altimetry  data,  which  would  be  inexpensive  to  implement, 
but  large  quantities  of  data  from  Navy  Navigation  Satellites  will  pro- 
bably have  to  be  processed  with  the  altimetry  data  to  achieve  improve- 
ment in  orbit  accuracies  currently  being  achieved.  Such  processing 
would  be  expensive.  A larger  number  of  gravity  coefficients  would  in- 
crease processing  costs  by  a low  or  at  worst  a moderate  amount.  Miscel- 
laneous changes  are  being  made  in  the  computations  to  account  for  the 
effects  of  atmospheric  and  ocean  tides  on  the  force  field  and  for  the 
effect  of  solid  earth  tides  on  station  position.  These  will  incur 
negligible  costs  in  operation  but  will  not  yield  significant  improve- 
ment until  gravity  and  instrument  errors  are  reduced.  Processing  data 
for  additional  satellites  will  increase  precision  by  an  amount  inversely 
proportional  to  the  square  root  of  a number  of  satellites  but  more  im- 
portantly will  shorten  the  time  span  required  for  a given  precision, 
thus  eliminating  sources  of  systematic  error  due  to  changes  in  station 
array.  During  some  time  periods,  these  calculations  have  been  made  but 
not  used  in  positioning.  Reprocessing  old  data  to  recompute  the  orbits 
would  reduce  a major  source  of  error  introduced  by  changes  in  coordinates 
of  base  stations  but  would  be  expensive.  Orbit  calculations  based  on 
drag-compensated  satellites  will  eliminate  another  source  of  error  in 
satellite  orbit  and  compute  pole  position. 


The  program  used  to  compute  station  positions,  "DARCUS,"  is 
undocumented  and  has  incurred  extensive  changes  over  the  original  code. 


Recoding  would  facilitate  configuration  control  and  give  greater  confi- 
dence in  the  results  but  would  be  expensive.  An  alternative  is  the  use 
of  another  program  that  is  available,  "CELEST  DARCUS,"  which  is  somewhat 
more  costly  in  operation  and  may  require  some  interfacing.  A tropospheric 
refraction  bias  parameter  exists  in  all  programs  at  NSWC,  and  its  use 
should  improve  heights  and  probably  longitudes  of  any  sites  for  which 
recomputations  are  performed  at  negligible  cost.  A correction  for 
second-order  ionospheric  refraction  could  improve  results.  Although 
simulations  made  to  date  have  not  indicated  that  the  correction  is 
significant  for  other  applications,  the  problem  should  be  reexamined 
for  this  purpose.  It  seems  likely  that  long  period  changes  in  solar 
activity  could  appear  as  long  period  changes  in  computed  station  height. 
The  first  order  ionospheric  refraction  recorded  by  the  geoceiver  would 
be  useful  for  this  purpose,  while  model  data  should  be  adequate  for  old 
data  and  other  base  stations.  Combining  observations  made  over  30-sec 
intervals  to  construct  60-sec  data,  after  first  filtering  the  30-sec  data 
to  detect  bad  data,  offers  hope  for  some  improvement  in  precision.  Going 
further  and  treating  the  data  as  biased  anges  in  the  manner  normally 
used  by  Brown6  yields  an  increase  in  precisions  of  a factor  of  three, 
which  would  be  very  valuable  if  the  oscillator  stability  is  sufficiently 
good  to  support  the  technique.  Both  of  these  approaches  are  being 
tested.  Consideration  of  the  orbit  error  in  the  solution  for  station 
coordinates  would  give  lower  elevation  passes  their  proper  weight  in  the 
solution,  although  common  orbit  biases  for  neighboring  stations  such 
as  employed  in  the  Brown6  short  arc  approach  appear  to  have  limited 
usefulness  for  the  generally  widely  separated  base  stations.  The 
comments  made  earlier  about  using  data  for  satellites  and  more  stations 
and  about  reprocessinq  old  data  apply  more  in  this  section,  except  that 
the  cost  is  small  compared  with  the  orbit  computation.  Some  recomputa- 
tions will  be  performed  in  any  case  to  isolate  the  cause  of  the  secular 
change  in  the  longitude  of  the  station  net. 

Survey  records  will  be  reexamined  in  a few  cases,  but  a continuing 
liaison  activity  would  be  required  to  insure  that  undocumented  changes 
are  not  made  to  antenna  installations  as  personnel  change  over  the 
relativley  long  period  of  time  that  would  be  required  for  Doppler 
determination  of  plate  motion.  Computations  of  rate  of  change  of  sta- 
tion position  could  be  improved  by  computing  relative  rates  of  change 
of  station  position  for  station  pairs  of  groups  over  identical  data 
spans.  This  decreases  the  precision  of  the  determination  but  increases 
the  accuracy.  The  difficulty  is  in  deciding  whether  to  consider  a 
surveyed  change  in  antenna  position  to  have  been  measured  perfectly  or 
to  be  a break  in  operation  of  the  station.  Even  if  realistic  error 
estimates  could  be  made  for  the  survey,  a solution  for  station  rate 
with  the  survey  connection  given  that  a-priori  uncertainty  would  probably 
result  in  a fit  to  errors  from  other  sources.  A limited  effort  will  be 
made  in  these  directions  as  part  of  the  evaluation  of  the  stability  of 
the  coordinate  system. 
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PLANS  AND  RECOMMENDATIONS 

A number  of  actions  are  being  taken  to  improve  the  accuracy  of 
determining  coordinates  of  stations  making  Doppler  observations  of  Navy 
Navigation  Satellites.  These  actions,  noted  by  "N"  in  the  fourth  column 
(Special  Effort  Required)  of  Table  5,  may  result  in  the  determination 
of  plate  tectonic  motion  from  these  observations.  To  increase  the 
probability  of  success,  obtain  results  for  more  geographical  areas,  and 
achieve  the  results  more  quickly,  the  actions  noted  by  "Y"  in  that 
column  of  the  table  should  be  taken. 

A decision  on  most  of  these  items  is  not  time-critical,  but  it  is 
impossible  to  totally  compensate  computationally  for  changes  in  the  base 
station  network.  Sites  such  as  Australia  and  Brazil,  which  are  distant 
from  other  base  stations  and  have  the  longest  period  of  operation  using 
the  same  antenna,  are  particularly  critical  to  validation  of  the  stability 
of  the  coordinate  system  and  of  great  importance  to  plate  tectonics 
studies.  Note  that  computed  latitude  rate  for  Australia  relative  to  the 
mean  coordinate  system  of  the  Doppler  base  stations  is  12.1+1.5  cm/yr 
over  the  last  5 yr  and  12.7+1.0  cm/yr  over  the  last  15  yr.  These  values 
are  consistent  with  the  motions  of  the  continent  reported  for  geologic 
time  periods.  It  is  recommended  that  every  reasonable  effort  be  made  to 
maintain  the  current  base  station  network  to  obtain  valid  data  for  other 
components  of  position,  other  geographic  areas,  and  other  time  periods. 
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ABSTRACT 


Doppler  observations  over  a three  year  period  on  polar  orbiting 
satellites  were  analyzed  to  determine  the  extent  to  which  the  movement 
of  tracking  sites  could  be  defined.  The  motion  of  sites  with  respect 
to  the  major  western  hemisphere  land  masses  could  only  be  determined 
to  an  accuracy  of  about  100  cm/year.  An  accuracy  of  5 cm/year  appears 
achievable  for  an  average  over  a 10-20  year  period  if  the  observation 
program  is  continued. 

BACKGROUND 

Runcorn's  studies  [1962,  1967a,  1967b]  have  intensified  interest 
in  relating  horizontal  motions  of  the  continents  to  large-scale  con- 
vection currents  in  the  earth's  mantle.  Supporters  of  the  relation- 
ship include  Allan  et  al  [1967],  Bott  [1967],  Orowan  [1967],  and 
Oxburgh  [1967].  Kaula  [1963]  and  MacDonald  [1964]  not  only  question 
the  validity  of  the  convection  hypothesis  but  they  and  other  scientists 
also  have  reservations  about  the  evidence  indicating  the  existence  of 
continental  drift,  which  is  required  to  support  the  convection  theory. 
Creer  [1967],  Girdler  [1965]  and  Bullard  et  al  [1965]  are  among  those 
who  have  advanced  data  supporting  continental  drift  while  Hospers  [1967], 
Lyustikh  [1967]  and  MacDonald  [1963]  have  proposed  alternate  inter- 
pretations of  the  data.  Satellite  observations  offer  a possible  addi- 
tional source  of  information  on  the  motions  of  the  continents.  Such 
data  have  been  reported  to  be  capable  of  positioning  the  observing 
station  to  10  m accuracy  [Anderle  1967,  Clancy  1967].  While  such 
accuracy  does  not  appear  to  be  adequate  to  resolve  the  question  of  the 
existence  of  a 1 cm/year  drift,  the  rapid  advancements  in  dynamic 
satellite  geodesy  have  made  it  appear  useful  to  evaluate  the  position- 
ing accuracy  of  the  most  recent  data  more  carefully. 

PROBLEM  AREAS 

This  study  is  based  upon  satellite  observations  made  by  the  U.  S. 

Navy  Doppler  system  which,  as  described  by  Newton  [1967],  has  an  in- 
strument accuracy  of  about  10  meters  for  two  dimensional  positioning 
of  the  observing  site  from  data  taken  on  a single  satellite  pass. 

Since  thousands  of  satellite  passes  are  observed  each  year  at  each 
instrument  site,  the  mean  station  position  could  be  considerably 
better  to  the  extent  permitted  by  computational  techniques.  The 
limiting  computational  factor  in  the  past  has  been  the  accuracy  of  the 
determination  of  the  gravitational  field  of  the  earth.  Uncertainties 
in  the  gravity  field  produce  errors  in  the  computed  satellite  orbit 
which  are  transferred  as  errors  in  the  computed  satellite  orbit  which 


are  transferred  as  errors  in  the  computed  position  of  the  observing 
site.  Additional  computational  errors  arise  from  the  neglect  of  the 
effects  of  solid  earth  tides  on  the  shape  of  the  earth  and  the  neglect 
of  the  effects  of  ocean  tides.  Effects  of  solid  earth  tides  on  the 
force  field  about  the  earth  are  only  included  in  the  most  recent  com- 
putations. The  principal  effect  of  these  computational  errors  will 
be  to  produce  random  error  and  a fixed  bias  in  relative  station  posi- 
tion. However,  should  relative  station  motions  be  found,  it  seems  un- 
likely that  conclusive  proof  could  be  could  be  established  that  the 
measured  motions  are  not  due  to  computational  errors.  On  the  other 
hand,  the  possibility  that  computational  errors  would  accidently  cancel 
any  actual  station  motions  seems  quite  remote.  Therefore,  analysis  of 
the  satellite  data  could  produce  evidence  refuting  the  existence  of 
continental  drift  at  the  present  time,  but  any  evidence  favoring  drift 
might  be  debatable.  An  additional  difficulty  in  the  use  of  satellite 
observations  for  the  intended  purpose  is  uncertainty  in  antenna  posi- 
tion. Since  two  antennas  separated  by  a few  meters  are  used  in  the 
correction  of  ionospheric  refraction  effects,  the  effective  station 
position  is  difficult  to  define  and  may  in  fact  vary  according  to  the 
geometry  of  a particular  pass  of  the  satellite  over  the  station. 

Errors  from  this  source  can  also  be  expected  to  affect  only  the  random 
error  of  observation  so  long  as  the  antennas  are  not  moved.  While  it 
is  possible  to  select  data  based  on  the  same  satellite  frequencies,  it 
is  not  a trivial  matter  to  have  assurance  that  the  precise  antenna  lo- 
cations are  unchanged.  In  the  early  development  of  the  Doppler  system 
antenna  experiments  might  have  caused  changes  in  station  position  ex- 
ceeding a meter.  Therefore,  the  older  observations  are  not  used. 

Even  currently,  if  construction  at  a station  requires  a change  in  an- 
tenna location,  the  accuracy  of  the  repositioning  may  not  exceed  30  cm. 
However,  the  record  of  antenna  positions  suggests  that  a repositioning 
accuracy  of  about  1 cm  can  be  achieved  for  several  sites. 

PROCEDURE 

The  source  data  for  these  calculations  were  normal  equations  for 
geodetic  and  orbit  parameters  formed  for  Doppler  Satellite  observations 
in  the  course  of  a general  geodetic  solution.  Each  normal  equation  was 
based  on  observations  made  in  either  an  8 or  16  day  time  span  at  13 
common  sites  and  at  up  to  5 additional  sites  which  were  usually  dif- 
ferent for  different  time  spans.  The  parameters  of  the  equations  in- 
cluded over  200  gravity  coefficients,  the  coordinates  of  all  the  ob- 
serving sites,  instrument  biases  (implicitly),  and  orbit  parameters 
which  included  the  six  orbit  constants,  one  for  solar  radiation  pres- 
sure, one  to  scale  the  drag  force  and  two  extraneous  parameters.  The 
extraneous  parameters  were  used  to  account  for  the  dominant  effect  of 
resonance  effects  due  to  gravity  coefficients  of  24th  to  28th  order. 

The  normal  equations  for  satellites  of  seven  orbital  inclinations  were 
combined  and  solved  to  produce  a set  of  geodetic  parameters  referred  to 
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as  NWL-8D.  For  purposes  of  this  study  the  individual  normal  equations 
for  polar  satellites  were  solved  separately,  employing  certain  con- 
straints, to  yield  solutions  for  8 to  16  day  time  spans.  The  solutions 
are  given  as  deviations  from  the  NWL-8D  solution.  Constraints  were 
generally  imposed  to  force  most  of  the  gravity  coefficients  and  the 
positions  of  stations  in  North  America,  South  America,  and  Greenland  to 
conform  with  the  NWL-8D  results.  The  second  degree  zonal  and  tesseral 
gravity  coefficients  were  generally  left  free  in  the  solution.  The 
effects  on  the  solution  of  small  changes  in  the  number  of  station  and 
gravity  parameters  were  tested  in  special  computer  runs  and  found  to  be 
negligible.  The  study  was  restricted  to  observations  of  polar  satel- 
lites in  order  to  use  homogeneous  data  for  the  maximum  available  time 
span.  While  data  for  nonpolar  satellites  were  available  for  earlier 
time  periods,  the  effects  of  remaining  errors  in  the  gravity  field  are 
different  for  satellites  having  different  orbital  inclinations  which 
would  lead  to  a spurious  difference  in  relative  station  positions  be- 
tween the  older  and  the  more  recent  observations. 

Corrections  for  changes  in  antenna  locations  were  made  on  the 
assumption  that  the  electrical  center  was  .163  of  the  distance  from 
the  400  me  antenna  to  the  150  me  antenna.  The  correction  is  based 
on  the  pair  of  frequency  equations 

f 1 = fTl  - — ^ + -f-U- 

C tT1 


f2  = f 2 “ f T2  + a°° 

C f T2 

where  fj  and  f2  are  the  frequencies  received  on  the  150  and  500  me 
antennas,  respectively. 

f , » 150  me, 

T 1 

f - 400  me, 
t2 

c - velocity  of  light, 

f and  fa,  are  the  range  rates  with  respect  to 
the  positions  of  the  150  and  400  me  antennas,  respectively,  and 
a^/^j  and  a^/f^  are  t*,e  fir8t  order  ionospheric  refraction  cor- 
rections at  the  positions  of  the  150  and  400  me  antennas,  respectively. 

An  analog  combination  of  these  two  frequencies  is  performed  at  the 
station  in  order  to  eliminate  the  refraction  term.  The  following  com- 


bination, or  a known  multiple  of  it,  is  recorded: 
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where  R » f /f  = 150/400. 

T 1 T 2 

The  first  two  terms  on  the  right  give  the  refraction  corrected  fre- 
quency expected  at  the  high  frequency  antenna  location  "B."  Setting 
+ jH?  (xy  -xj  , the  third  term  yields  the  correction 
3x 

R2/ (1  - R2)  of  the  distance  from  location  B to  location  A.  The  final 
term,  the  difference  in  refraction  at  the  two  positions,  is  neglected. 

RESULTS 

A total  of  25  separate  solutions  for  station  coordinates  was 
computed  using  polar  satellite  observations  made  in  8 or  16  day  time 
spans  between  late  1963  and  early  1967.  The  orbital  elements  for  the 
different  satellites  are  given  in  table  1.  Figure  1 shows  the  solu- 
tions obtained  for  the  site  in  Misawa,  Japan,  holding  the  positions  of 
sites  in  Greenland,  North  America  and  South  America  fixed  at  values 
found  in  the  NWL-8D  solution.  A trend  in  the  longitude  of  the  station 
can  be  noted,  although  the  statistical  significance  of  the  trend  is  not 
apparent.  In  order  to  provide  guidance  on  the  statistical  significance 
of  the  results,  a least  squares  fit  of  a straight  line  was  made  to  the 
solutions  for  each  coordinate  of  each  permanent  station.  Solutions 
based  on  fewer  than  21  satellite  passes  were  ignored.  The  rms  resid- 
uals, the  slope  and  the  standard  error  in  the  slope  are  given  in 
table  2 and  graphed  in  figure  2.  The  drift  of  the  stations  is  gen- 
erally not  significant  compared  to  the  accuracy  of  the  determination 
which  is  about  100  cm/year  for  the  stations  which  made  the  most  com- 
plete set  of  observations.  Although  there  are  a few  cases  where  the 
computed  drifts  exceed  the  standard  error  in  the  drift,  the  biases  are 
probably  due  to  a defect  in  the  statistical  model:  The  individual  so- 
lutions were  assumed  to  be  of  equal  quality  while  in  fact  the  older 
solutions  appeared  to  have  a larger  scatter  than  those  based  on  more 
recent  and  more  numerous  solutions.  The  overall  standard  deviation  of 
an  individual  solution  of  about  4 meters  is  a factor  of  two  or  three 
better  than  that  reported  a few  years  ago  [Anderle  1967]. 

SUMMARY 

The  accuracy  obtained  in  the  determination  of  station  movement  is 
not  sufficient  to  either  support  or  refute  the  conclusion  that  the 
continents  are  drifting.  However,  the  accuracy  achieved  is  better 
than  has  generally  supposed  to  be  possible.  Since  the  accuracy  of  the 
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drift  determination  varies  linearly  with  the  time  span  of  observation, 
an  increase  in  time  span  of  observation  from  the  present  3 years  to  20 
years  would  improve  the  accuracy  of  the  determination  to  15  cm/year. 

The  present  study  is  based  on  about  six  months  of  data  observed  on  one 
or  another  satellite  during  the  three  year  period  considered;  if  the 
same  proportion  of  data  were  processed  for  the  20  year  period,  the 
accuracy  would  be  further  improved  to  about  5 cm/year.  Still  better 
accuracy  could  be  achieved  by  virtue  of  a modest  additional  improvement 
in  accuracy  or  by  the  use  of  a larger  proportion  of  the  data  obtained. 
Therefore,  meaningful  results  would  probably  be  obtained  in  about  20 
years  if  the  observation  program  continues. 
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TABLE  1 

PRINCIPAL  ORBITAL  ELEMENTS 


Satellite 

Inclination 

Perigee  Height (km) 

Eccentricity 

1963  38C 

89.9° 

1067 

.004 

1963  49B 

90.0° 

1056 

.003 

1964  26A 

90.5° 

864 

.007 

1965  48A 

90.0° 

1025 

.007 

1965  109A 

89.1° 

912 

.010 

1966  24A 

89.7° 

888 

.016 

1966  41A 

89.9° 

856 

.008 

1966  76A 

88.8° 

1048 

.003 
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Figure  2h.  Station  Coordinate  Changes  (Greenland) 
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Figure  2m.  Station  Coordinate  Changes  (England) 
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Figure  2s.  Station  Coordinate  Changes  (Phillipines) 
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Figure  2u.  Station  Coordinate  Changes  (Samoa) 


suim 


n Kiumi 


( »u»») 


LOMI  T U 0 1 < L ) 

l*i  i mod  f i 

RUBRIC  R> 

rnttei  » 

milt  o 

Mii 

IS 

11. B 

1 1 .0 

11.4 

mums 

•crm 

Ml 

CM 

10?  .0 

• 1.0 

!»!.♦ 

• MM 

t R/OB 

•1M 

1 .1 

-II  .1 

IMOI 

CR/OB 

M 

0 . 1 

4.4 

BCllBRBlS 

• MKB 

Ml 

CR 

01. 1 

41. 1 

OB. 4 

N «» 

- t H 


JjJ 

V 


®w 

/V^V-,  4*- 


■V 

■"  .an  1. 


Figure  2x.  Station  Coordinate  Changes  (McMurdo) 
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Station 

Longitude 

Latitude 

Height 

North  America 
California 

New  Mexico 

6.03.5 

15.1+3.0 

-7 .5+2.0 
-8. 3+2. 2 

16.6+2.2 
-1.3+3. 3 

South  America 
Brazil 

-5. 4+3. 7 

-2. 6+3.0 

-15-1+3.8 

Europe 

England 

Belgium 

-8. 7+2. 8 
-25.9+2.6 

0. 9+2.0 

8. 3+2.0 

34.6+3.0 

20.6+2.5 

Africa 

South  Africa 
Seychelles^1' 

-7.1+3. 3 
-17-5+4.0 

-1.4+2. 2 

3. 5+2.6 

-18.4+3.3 

30.4+3.0 

Asia 

Japan 

Phillipines 

-19-4+3.2 

-15.4+4.2 

-11.4+3.1 

-18.9+3.5 

25.3+3.4 

10.3+3.0 

Pacific 

Hawai ( 1 ) 

-3. 3+4. 2 

-1.8+2. 6 

7 . 4_3 . 0 

Australia 

Adelaide 

-27.1+2.1+ 

12.1+1.5 

2. 7+1. 7 

Antarctica 

McMurdo 

-33.9+5-7 

-1.2+6. 3 

7. 2+5. 8 

RMS 

12. 5^2^ .3.6 

8.4, 3.0 

18.8,3.6 

(1) Solutions  were  not  computed  in  1973. 

( 2 )  About  mean 


Table  2a.  Position  Change  Rates  (Longitude,  cm/yr) 
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H.A.  - Not  applicable  or  not  available  due  to  short  data  span. 

(1)  - Station  va3  not  actually  moved  but  different  antennas  arc  used. 

(2)  - Sane  as  1972.0-1978.0  because  station  started  operations  after  1972.0. 


R.A.  - Rot  applicable  or  not  available  due  to  short  data  span. 

(1)  - Station  vas  not  actually  acrved  but  different  antennas  are  used. 

(2)  - Sa*e  as  1972.0-1978.0  because  station  started  operations  after  1972.0. 


Table  2c.  Position  Change  Rates  (Height,  cm/yr) 
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N.A.  - Sot  applicable  or  not  available  due  to  short  data  span. 

(1)  - Station  vas  not  actually  moved  but  different  antennas  are  used. 

(2)  - Same  as  1972.0-1970.0  because  station  started  operations  after  1972.0 
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Table  3.  Station  Coordinate  Changes 
Over  14  Years  (cm/vr^ 

STATION  LONGITUDE  LATITUDE  HEIGHT  NUMBER  OF  SOLUTIONS 


NORTH  AMERICA 


Maine!^)  . . 

6. U+l. 8 

15-9+1.5 

2. 9+1- 6 

331 

Minnesota^ 2 ' 

-22.8+2.0 

-5. 9+1. 3 

-lU.8+1.9 

398 

California 

2. 5+1. 6 

-1. 3+1.1 

1*. 7+1.0 

373 

Maryland! » (2) 

U.U+2.0 

U . U+l . u 

lU.0+2.7 

259 

Alaska!1) • (2) 

10.1+1.6 

-9. 0+1.2 

29.6+1.8 

237 

Greenland ( 1 ) » ( 2 ) 

-10.5+1.1 

6. 2+1.1 

19.2+1.6 

296 

SOUTH  AMERICA 

Brazil 

-3.U+2.1 

-7. 6+2.0 

-9. 9+2. 2 

3U7 

AFRICA 

South  Africa 

2. U+l. 8 

-10.6+1.3 

10.5+1.9 

378 

-8. 6+2.1 

0.7+0. 6 

-6. U+l. 8 

322 

ASIA 

Japan 

1.9+1. 8 

0.1+1. 6 

21.3+1.7 

UlO 

Phillipines 

-6. 3+2.0 

-7. 2+2.0 

9. 2+2. 5 

313 

PACIFIC 

Hawaii 

3. 8+1. 9 

9. 8+1. 2 

-5. 9+1- 2 

328 

AUSTRALIA 

Adelaide 

-lU.2+1.3 

12.7+1.0 

21. U+l. 0 

397 

ANTARCTICA 

McMurdo  -23.713.1*  -0.2+3. 8 60.1+3.5  286 

RMS  11.0,1.9  8. 1,1. 7 15.1(3),1.6 


NOTES : (l)  1973  data  not  processed 

(2)  station  closed  or  moved  one  to  two  years  before  end  of  data  span 

(3)  excluding  Antarctica 
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Table  4.  RMS  Residuals  from  NWL  9Z-2  Coordinates  (cm) 


NORTH  AMERICA 

Longitude 

Latitude 

Height 

Maine 

no 

106 

1U8 

Minnesota 

130 

63 

106 

California 

93 

52 

Sh 

New  Mexico 

87 

6i 

93 

Maryland 

99 

57 

132 

Texas 

118 

63 

121 

Alaska 

92 

U5 

7U 

Greenland 

66 

56 

88 

Ottawa 

119 

73 

lkl 

Calgary 

121 

118 

96 

Virginia 

90 

137 

86 

SOUTH  AMERICA 


Brazil 

99 

77 

115 

EUROPE 

England 

97 

56 

197 

Belgium 

8U 

57 

79 

Italy 

97 

66 

6>4 

AFRICA 

South  Africa 

101 

65 

105 

Seychelles 

107 

65 

97 

ASIA 

Japan 

100 

93 

109 

Phillipines 

118 

97 

151 

PACIFIC 

Hawaii 

108 

67 

86 

Samoa 

llU 

102 

176 

Guam 

107 

6U 

193 

AUSTRALIA 

Adelaide 

96 

<*9 

52 

ANTARCTIC 

McMurdo 

129 

136 

20*4 

AVERAGE 

103 

76 

116 

C-6 
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Table  5.  Potential  for  Improvement 


3.  Orbit 

Gravity  field 

Misc  orbit  changes  (ocean,  atmrsph., 
land  tides) 

Process  more  satellites 
Re-process  old  data 
Use  drag-free  satellites 

1.  Rubidium  Oscillator  for  Geoceiver 
Reprogram 

Tropo  bias  parameter 

2nd  order  ionospheric  correction 

Longer  count  interval 

Biased  range 

Orbit  bias 

Process  more  satellites 
Process  available  data  for  more 
stations 

Re-process  old  data  - New  orbits 

- Current  orbits 

5.  Review  Survey  Records 

6.  Rate  Computations 

Simultaneous  treatment  of  data 
Consider  survey  connection  uncertainty 

2.  Data  Acquisition 

Encourage  more  stations 
Maintain  fixed  network 


M 


N 


N 


M 

H 

M/N 

M 

M 

L 

M 


H 

N 

M/L 

L 

L 

M 

N 


N 

N 

M/L 

L 

L 

L 

M 


N 

N 

Y 

Y 

Y 

Y 
N 


Y 

I 

Y 

Y 

Y 

Y 

Y 


M 

M 

M 


Y 

Y 

Y/N 


Y/N  N 


H 

M 


Y/N  N 
Y/N  N 


L 

H/L 


N 

N 


Key:  H,  M,  L are  high,  medium,  low. 
Y,  N are  yes  and  no  or  none. 
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DISTRIBUTION  (cont'd) 


HQ  SAMSO 

Air  Force  Unit  Post  Office 
Los  Angeles,  CA  90045 
Attn : YE 

U.  S.  Naval  Academy 
Annapolis,  MD  21402 
Attn:  Technical  Library 

Chief  of  Naval  Research 
Department  of  the  Navy 
Arlington,  VA  22217 
Attn:  Code  418 

Chief  of  Naval  Operations 
Department  of  the  Navy 
Washington,  DC  20350 
Attn : OP986 

Director 

U.  S.  Army  Ballistic  Research  Laboratory 
Aberdeen  Proving  Ground,  MD  21005 
Attn:  DRDAR-TSB-S  (STINFO) 

Local: 

E41 
K10 
Kll 
K12 
K13 
X210 

X2101  (GIDEP) 


(100) 

(10) 

(10) 

(10) 

(2) 


